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ABSTRACT 

We present 870 observations of dust continuum emission from the LkCa 15 pro- 
toplanetary disk at high angular resolution (with a characteristic scale of 0'.'25 = 35 AU), 
obtained with the IRAM Plateau de Bure interferometer and supplemented by slightly 
lower resolution observations from the Submillimeter Array. We fit these data with sim- 
ple morphological models to characterize the spectacular ring-like emission structure of 
this disk. Our analysis indicates that a small amount of 870 /xm dust emission (~5 mJy) 
originates inside a large (40-50 AU radius) low optical depth cavity. This result can be 
interpreted either in the context of an abrupt decrease by a factor of ~5 in the radial 
distribution of millimeter-sized dust grains or as indirect evidence for a gap in the disk, 
in agreement with previous inferences from the unresolved infrared spectrum and scat- 
tered light images. A preliminary model focused on the latter possibility suggests the 
presence of a low-mass (planetary) companion, having properties commensurate with 
those inferred from the recent discovery of LkCa 15b. 

Subject headings: circumstellar matter — protoplanetary disks — planet-disk interac- 
tions — submillimeter: planetary systems — stars: individual (LkCa 15) 



1. Introduction 

Hundreds of exoplanets have been discovered around main-sequence stars, a nd substantial ef- 



fort is being invested to explain their demographics with formation models (e.g.. Ilda &: Linll2004l : 



Mordasini et al.ll2009l ). But associating exoplanet properties with their formation epoch is prob- 
lematic: dramatic evolutionary processes that occur at early times are closely tied to the unknown 
physical conditions in the progenitor circumstellar disk. Ideally, mature exoplanets could be com- 
pared with their younger counterparts, still embedded in their natal disks. However, detecting 
planets around young stars is difficu lt. Radial velocity and transit searches are hindered by stellar 



variability (e.g.. lHuelamo et al.ll2008l ). and direct imaging is limited by contrast with the bright star 
and disk emission. However, the presence of a young planet can be inferred indirectly through its 
dynamical imprint on the structure of the disk material. A sufficiently massive planet (> 1 Mj^p) 
opens a gap that impedes the inward flow of mass through the disk, decreasing the densities at the 
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disk center (e.g.. lLin &: Papaloizoul 119861 : iBrvden et al.l 11993 : iQuillen et al.ll2004l ). The location of 
the gap mar ks the planet orbit, and the amount of materi al that flows across it depends on the 
planet mass (jLubow &: D'Angeldl2006l : IVarniere et alj 120061 ) . In principle, the orbit and mass of a 
~Myr-old giant planet can be estimated from observations of its disk birthsite, through constraints 
on the gap location and the amount of material interior to it, respectively. 

The disk around the young star LkCa 15 is considered an excellent candidate for planet- 
induced disk clearing, based on its distinctive infrared spectrum (lEspaillat et alJ 120071) as w ell as 
the ring-like morphology of its mm- wave dust einission ( Pietu et alj|2006l : [Andrews et al.ll201ll ) and 
scattered light in the infrared (jThalmann et al.ll20ld ). Those observations confirm that the LkCa 
15 disk has a large central "cavity", with significantly diminished dust optical depths on Solar 
System size-scales. However, the cavity is not empty. A faint infrared signal is detected in excess of 
the stellar ph otosphere, indicating that at least a small amount of warm dust resides near the star 
(jEspaillat et al. 2008). That excess verifies the presence of a tenuous inner disk - and therefore a 
gap - although it provides only minimal bounds on its size (and therefore the gap width) and mass. 
Based on an attempt to model a h igh resolution Submillimeter Array (SMA) observation of the 
LkCa 15 disk. [Andrews et al.l (j201ll ) identified preliminary evidence for weak, optically thin 870 /;xm 
emission from dust inside the disk cavity. If confirmed, that emission can be used to estimate the 
inner disk mass, a key diagnostic of the fiow rate across the gap. 

In this Letter, we present new 870 /xm continuum observations of the LkCa 15 protoplanetary 
disk, with a 50% improvement in angular resolution facilitated by the recent commissioning of 
high-frequency receivers at the Plateau de Bure interferometer (PdBI). In §2, we provide a brief 
overview of the new data and describe how their combination with previous SMA observations 
provide the sharpest view yet of the thermal emission from the LkCa 15 disk. In §3 we use simple 
models to explore the properties of the disk cavity and its contents. And in §4 we discuss those 
modeling results in the contexts of planet formation around LkCa 15 and the potential future utility 
of similar observations as an independent check on the properties of young exoplanets. 



2. Observations and Data Reduction 

LkCa 15 was observed for 5 hours with the most extended configuration (A: baselines of 
130-760 m) of the PdBI on 2011 January 27. The observations were conducted in "shared-risk" 
mode since they used the new Band 4 receivers at an effective continuum frequency of 345.8 GHz 
(868 fxm) and the new WideX correlator to sample the continuum emission with a total bandwidth 
of 3.6 GHz (per polarization). The observations cycled between LkCa 15 and two nearby quasars, 
J0530+1331 and J0336+3218, every 22 minutes. The data were calibrated with the CLIC software 
in the GILDAS package. Short observations of the bright quasars 3C 454.3 and 3C 273 were used 
to set the bandpass and absolute flux scale, and the nearby quasars that were interleaved in the 
observing cycle were utilized to calibrate the time-dependent complex gain response of the system. 
At the time of the observations, the new Band 4 LO system perturbed the first channel (of 3) in 
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the PdBI water vapor radiometer (WVR) phase correction system. We reduced the WVR system 
to a dual channel mode in the post-processing, and smoothed the WVR data on 5 s intervals. The 
differential phase correction determined on 45 s intervals was extended over each source cycle by 
fitting and removing linear instrumental drifts. This process requires a stable atmosphere, with 
water vapor fluctuations that average to near zero over the source cycle. These conditions were 
generally met, due to the low water vapor levels (<2mm) present throughout the observations. 

To improve the Fourier cover age on short spacings , we supplemented these PdBI observations 
with the SMA data described bv [Andrews et al.l (|201ll . baseline s of 8-508 m). After adj usting the 
datasets to account for the small proper motion of LkCa 15 (jPucourant et al.ll2005l ). the disk 
cent roid was estimated i n each dataset by minimizing the imaginary components of the visibilities 
(see lAndrews et al.ll201ll ). The inferred reference centers for the two datasets agree within ~10mas 
and are <70mas from the expected stellar position (within the absolute astrometric uncertainty in 
each dataset), at RA = 4'^39"^ir;80 and DEC = -h22°21'03'.'20. The SMA and PdBI calibrations 
were compared over their redundant Fourier coverage, and were found to be in excellent agreement 
on 150-500 kA baselines: deviations between the visibility amplitudes in each dataset are random, 
with an RMS difference of <5%. The combined SMA and PdBI visibilities were Fourier inverted 
assuming natural weighting, deconvolved with the CLEAN algorithm, and restored with a 0'.'33 x 0'.'22 
synthesized beam using the MIRIAD software package. The resulting synthesized continuum map 
is shown in Figure [H with an effective wavelength of 870 //m, RMS noise of 0.7 mJy beam~^, peak 
flux density of 27mJy beam~^, and integrated flux density of 380 mJy. 



3. Results 



The 870 ^um image in Figure [T] provides the shar pest view yet of cool dust emission fro m the 
LkCa 15 disk. As noted previously at lower resolution ( Pietu et al. 2006 : Andrews et al.l[201ll ). this 



emission has an inclined ring morphology with a large and prominent central depression in intensity. 
The emission ring peaks at semimajor separations of '^0'.'4 (56 AU for an assumed distance of 140 pc) 
and has an aspect ratio and orientation in good agreement with the inclinat ion (i = 51°) an d 



major axis position angle (PA = 61°) inferred from its molecular line emission (jPietu et al.l 120071 ). 
Figure [2] shows the azimuthally- aver aged visibilities as a function of the deprojected baseline length 
(accounting for the disk viewing geometry) . The real part of this visibility profile exhibits the classic 
oscillation pattern expected from the Fourier transform of a ring in the sky-plane, with distinct 
nulls (sign changes) at deprojected baselines near 150, 350, and 700 kA. The imaginary terms 
are negligible on all baselines, consistent with an axisymmetric emission distribution. Although 
subtle, two qualitative features in the data can serve as useful benchmarks in a refined effort to 
characterize the LkCa 15 disk structure. First, the continuum intensities inside the ring are small, 
but not zero (see Figure [1]). And second, the oscillations in the continuum visibility profile are 
relatively muted, with a maximum amplitude of only ~5 mJy between the second and third nulls. 
This latter property suggests that the emission peak near the inner ring edge is not very sharp. 
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With those features in mind, we attempted to reproduce these LkCa 15 disk observations with 
simple 870 /um emission models. We adopted a radial surface brightness prescription that assumes 
optically thin thermal emission, ly oc B^{Td){l — e""^) ~ By{Tci)T, where Bj^ is the Planck function, 
T(i the dust temperature, and r the optical depth. The temperature profile was fixed to Td{R) = 
100(i?/l AU)~'^'^ K, based on a crude approxir aation of the midplan e temperatures derived in a more 



sophisticated treatment of radiative transfer (jAndrews et al.ll201ll ). Assuming the dust emissivity 



is independent of radius, we utilized a parametric form for the base optical depth profile motivated 
by th e surface densities in idealized viscous accretion disks: rfc(i?) oc {R/ Rc)'"^ exp[—{R/ Rc)'^~'^] 



(e.g., iHartmann et al.lll998l ). Modifications to that base model were also considered, including an 
optical depth cavity where t{R < -Rcav) = ^Tb- Three model permutations were investigated: the 
base model {6 = 1, Rcav undefined; Model A), the base model with an empty cavity {6 = 0; Model 
B), and the base model with a partially filled cavity (0 < 6 < 1; Model C). All models have three 
base parameters - a gradient (7), characteristic size (Re), and normalization (defined as the flux 
density, Ftot = / lydVL) ~ and can utilize up to two additional parameters, {i?cav) <^}- 

For a given model type and parameter set, synthetic visibilities were computed for the appro- 
priate viewing geometry at the spatial frequencies observed by the SMA and PdBI. Those model 
visibilities were compared with the data and assigned a fit quality statistic, the sum of the (real and 
imaginary) values over all spatial frequencies. The best-fit parameter values for a given model 
were determined by minimizing with the Metropolis algorithn i, utilizirig mul tiple Monte Carlo 
Markov chains and an initial period of simulated annealing (see iGregorvl l2005l ) . The results are 
compiled in Tabled! The estimated parameter uncertainties do not consider correlated errors from 
the (fixed) temperature profile or viewing geometry, and therefore are clearly under-estimated. The 
best-fit synthetic data products for each model type are directly compared with the observations 
in Figures [2] and [3l The corresponding radial brightness profiles are shown together in Figure HI 

For Model A, the observed emission morphology can only be rep roduced with a large and 



negative optical depth gradient parameter, 7 (e.g., see llsella et al.ll2009l ). The Model A fit does a 
relatively poor job accounting for the breadth of the observed ring structure: there is a tendency 
to over-predict the emission in the disk center and prematurely cut off at larger radii. Significant 
improvement is made with Mo del B, when a cayity is added to the base model. This is effectively 



the same structure assumed by lAndrews et al.l (|201ll ). That preliminary work used a fixed 7 = 1 



which tends to maximize the peak-to-cavity emission contrast in the fits, leading to higher positive 
residuals at the disk center. Similar results were obtained when that effort was repeated here, with 
strong centralized residuals (~11 a = 8.2 mJy). The best-fit Model B parameters are different than 
the fixed-gradient case (see Tabled! x\-xl = 1350 - but those same residuals remain significant 
(~7cT = 4.8 mJy). Naturally, this motivated the addition of an emission component inside the disk 
cavity, cast for simplicity as an adjustment to 5 (Model C). The inclusion of that weak emission 



^The differences in our progression of models are large enough (the best-fit likelihood ratios are significantly 
greater than unity) to warrant the complexity of adding a parameter at each step from Models A through C. 
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improved the fit quality {xb ~ Xq ~ leaving no significant residuals compared to the data. 

In this scenario, dust inside the disk cavity produces ~5 mJy of 870 /im emission, corresponding to 
20% of the peak surface brightness and only 1% of the integrated flux density. 

A gap structure represents an alternative model that naturally produces dust emission inside a 
disk cavity. To explore that possibility with a more ph ysically motivated p resc ription, we modeled 
the da ta with the treatment of gap profiles advocated bv lCrida et al.l (j2006l ) and lCrida &: MorbideUi 
(|2007l ). In this scenario (Model D), we utilized a semi-analytic approximation for the surface 
density perturbation produced by an embedded low-mass companion to modify the base optical 
de pth profile. The d epth, width, and shape of the gap profile perturbation were characterized 
bv ICrida et al.l (|2006l . their Eq. 14) in terms of the companion-to-star mass ratio {q = Mg/M^), 
the semimajor axis of the companion {Rg), the disk viscosity (z^), and the local disk aspect ratio 
{H/R, where H is the vertical scale height of the gas). Following Cr ida and his colleag ues, we fixed 
H/R = 0.05 and only investigated models where u = in the ICrida et al.l ( 2000 ) normalized 
units (for our fixed Td profile, this corres ponds to a typical viscosity coefficient a ~ 0.001 in the 
formulation of lShakurafc SunyaevI Il973). Furthermore, we fi xed Rg = 16 AU, in line with the 

see §4). With these simplifying 
The Model D structure also has 
116, comparable to Model C). 



2011 



recent detection of a faint companion (|Kraus &: Ireland 
assumptions. Model D has four parameters, {7, Rc, -Ftot, q}- 
improved fit quality relative to the empty cavity model (xb 
The estimate of q implies a companion mass of Mg = 9 it 1 M.Tnpi given t he Lk Ca 15 stellar mass 
of AI^ = 1.01 lb 0.03 Mq that was determined dynamically bv lPietu et al.l (|2007l ). We should again 
caution that these represent formal parameter uncertainty estimates that are only applicable under 
the restrained assumptions of this particular model: the true uncertainties could be significantly 
larger. As for Model C, there is roughly 5 mJy of 870 ^um emission interior to the gap of the favored 
Model D structure. 



4. Discussion 



We have used high angular resolution 870 fim PdBI+SMA observations to investigate the radial 
distribution of cool dust in the LkCa 15 protoplanetary disk with simple emission models. Although 
grounded in more sophisticated techniques, these models are inherently more morphological than 
physical. Their advantage lies in computation speed, which facilitated a broader exploration of dust 
structures that would have been prohibitive for a complex radiative transfer analysis. Despite their 
limitations, these simple models provide some fundamental qualitative insights on the LkCa 15 
disk properties: (1) there is a substantial decrease in the dust optical depths inside R ^ 40-50 AU; 
(2) the emission just outside that cavity edge is not sharply peaked, as attested by the smooth 
intensity profiles produced by the favored negative optical depth gradients (7); and (3) there is 
a small amount of dust located inside the disk cavity. Given our limited resolution, the spatial 
distribution of that weak emission in the cavity is unclear. It may fill the cavity (Model C), or it 
may be more centrally concentrated in the form of a gap structure (Model D) similar to what was 
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inferred from models of the unresolved infrared spectrum ([Espaillat et al.l 120081 ) 



If the latter is true, the gap is most likely opened by the resonant torques generated by 



1986 



Brvden et al 



inter actions between the disk and a low-mass companion (jLin &: Papaloizou 
19991 ) ■ Alternative gap-opening mechani sms - for example, photoevaporation - are unlik ely given 
the properties of the LkCa 15 system (jAlexander &: Armitagd l2009l : iQwen et alJ l201ll ). High- 
contrast imaging has ruled out stellar and brown dwar f companions around LkCa 15, hinting 
that the gap may be opene d by a young giant pla net (jThalmann et al.l l20ld : iPott et alJ l20ld : 
Kraus et alJ l201ll ) . Recently, iKraus &: Ireland! (j201ll ) used a non-redundant masking technique to 
detect a faint, co-moving companion ~0'.'07 from LkCa 15. If that object is co-planar with the disk 
and on a circular o rbit, it has a semimajor axis of 16 AU. Using a simple emission model based on 
the prescription of ICrida et alJ (|2006l ). we have shown that a gap at this location can reproduce 



well the resolved 870 Mm emission r norph ology we observe if the companion mass is ~9Mjup. At 



ages of 1-3 Myr, the iBaraffe et al.l ([2003) evolution models suggest that this object should have 
an infrared contrast of AX = 6.4-7.2, i n reasonable a greein ent with the IS.K = 6.8 measured by 
Kraus Ireland! (j201ll ). However, the iMarley et al.l (|2007l ) models suggest it would be ~150x 



fainter: a substantial accretion luminosity would be required to account for the observed infrared 
emission. 

Ultimately, improved constraints on the companion mass could be based on the disk contents 
interior to the gap. A crude estimate of the dust mass in that region can be made from the 
luminosity of the optically thin 870 ^um emission that was inferred in Models C and D. Assuming a 
dust opacity of 3cm^ g~^ and a fiducial Td = 45 K, the estimated flux density of 5mJy corresponds 
to 10~^ Mq (0.4 Mq). If that dust traces the gas at a rnass fr action of ~1%, then the accretion rate 
onto LkCa 15 (M* 2 x 10^^ vr~^: llnglebv et al.ll2009l ) implies that this inner disk material 
would rapidly drain onto the star (in <0.05Myr). Given the system age of 1-3 Myr, the inner 
disk must be continually replenished from the massive reservoir outside the gap. There is some 
notable tension with theoretical expectations here: it is not clear how a ~9Mjup companion can 
be reconciled with the inferred inner disk mass and stellar accretion rate in numerical simulations 



of gap-crossing flows (jLubow et al.l Il999l : iLubow &: D'Angeld 120061 ). If LkCa 15b has a much 
lower mass, it likely cannot sculpt the deep, wide gap needed to explain the obser v ation s: an 
additional companion with a lon ger orbital period must also be present. IZhu et al.l (120111 ) and 
Dodson-Robinson k. SalykI (j201ll ) have effectively argued for this latter possibility. They suggested 
that multi-planet systems can alleviate the apparent discrepancy between large transition disk 
cavities and accretion rates, implying that LkCa 15b is but one component in a young planetary 
system. 

Robust, quantitative constraints on the properties of LkCa 15b based on the structure of the 
LkCa 15 disk requires more work, including a stronger link between numerical simulations, an 
improved modeling effort, and observations that can probe the inner disk at even higher angular 
resolution. Nevertheless, the PdBI+SMA data presented here offer a tantalizing foreshadowing of 
the new roles mm-wave observations of disk structures can play in exoplanet science. 
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Table 1. Model Parameters 



Model ABC D 

Ftot (mJy) 363 ± 2 373 ± 2 367 ± 3 385 ± 2 

7 -1.7±0.1 -l.OiO.l -0.5±0.1 -0.3±0.1 

i?c (AU) 107 ±2 113 ±1 114 ±1 113 ±1 

i?cav (mJy) • • • 36 ± 1 49 ± 1 

6 1 (fixed) (fixed) 0.18 ±0.02 

R, (AU) • • • • ■ • • • ■ 16 (fixed) 

q ■■■ ■■■ ■■■ 0.009 ±0.001 

"x^ 516,735 516,600 516,489 516,484 



Note. — Parameter estimates, formal uncertainties, and values 
for the models discussed in §3. There are 776,966 independent visibil- 
ities used in the model fits. 
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Fig. 1. — Aperture synthesis image of the 870 fim continuum emission from the LkCa 15 disk, made 
from the naturally-weighted combination of PdBI and SMA datasets. The synthesized beam, with 
dimensions of Cf.'SS x [f.'22 (46 x 31 AU), is shown in the lower left. The wedge on the right marks 
the conversion from color to surface brightness. Each side of the image corresponds to 560 AU 
projected on the sky. 
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Fig. 2. — The real and imaginary 870 /xm visibilities as a function of baseline length, deprojected 
to account for the LkCa 15 disk viewing geometry and azimuthally averaged. The inset in the 
top panel is a detailed view of the gray-filled region. The best-fit models visibilities for different 
emission prescriptions are overlaid in color (all models have zero imaginary fluxes, by definition). 
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Fig. 3. — Comparison of the data and models in the image plane. The top left panel shows the 
same image as in Figured! To the right, the top panels display the best-fit model images, and the 
bottom panels the imaged residual visibilities. All panels show the same color scale and contour 
levels, starting at 1.4 mJy beam~^ (2<7) and increasing in 2.5 mJy beam"^ (3.5 cr) increments. As 
noted in Figure [21 Models C and D - which emulate a low-density (but not empty) cavity and a 
gap structure for the LkCa 15 disk, respectively - provide the best matches to the data. 




Fig. 4. — Radial surface brightness profiles for the best-fit parameters of each model type, cast 
for simplicity into a brightness temperature format. The combined PdBI+SMA data provide a 
maximum projected radial resolution of ~17AU, marked here by the shaded gray region. 



